Zero-field Kondo splitting and quantum-critical transition in double quantum dots 
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Double quantum dots offer unique possibilities for the study of many-body correlations. A system 
containing one Kondo dot and one effectively noninteracting dot maps onto a single-impurity An- 
derson model with a structured (nonconstant) density of states. Numerical renormalization-group 
calculations show that while band filtering through the resonant dot splits the Kondo resonance, the 
singlet ground state is robust. The system can also be continuously tuned to create a pseudogapped 
density of states and access a quantum critical point separating Kondo and non-Kondo phases. 

PACS numbers: 72.15.Qm,73.63.Kv,73.23.-b 



i 

C 

o 
o 



> 

00 

o 

oo 

o 
o 

-I— > 



i 

C 

O 

o 



The Kondo effect, arising from antiferromagnetic cor- 
relations between an unpaired spin and an electron bath 
Q , can be strongly modified by structure in the host den- 
sity of states (DoS). Geometric confinement and nar- 
row bands can dramatically change the Kondo state, 
with important observable consequences Q. In pseudo- 
gapped hosts, where the DoS vanishes as a power-law at 
the Fermi energy, a quantum critical point (QCP) sepa- 
rates the Kondo phase from one at smaller couplings in 
which the Kondo effect is completely suppressed [3, |y, Lj . 

Semiconductor quantum dots provide many opportuni- 
ties for systematic investigation of strong-correlation ef- 
fects . Single quantum dots have allowed controlled re- 
alizations of the Kondo regime of the Anderson impurity 
problem [jj. Recent attention has focused on the fasci- 
nating phvsics promised by double quantum-dot (DQD) 
systems For example, DQD experiments have inves- 
tigated the effect of interdot "hybridization" on Kondo 
physics ^lj , and have beautifully demonstrated the com- 
petition between the Kondo effect and the Ruderman- 
Kittel-Kasuya-Yosida interaction among localized spins 
|l2| . DQD setups have also been proposed to realize the 
unusual non-Fermi-liquid properties associated with the 
two-channel Kondo effect |l3). 

In this paper, we propose DQDs as a versatile exper- 
imental realization of an impurity coupled to an elec- 
tron bath having a structured (nonconstant) DoS. De- 
vices with one dot ("dot 1") in the Kondo regime and 
the other ( "dot 2" ) close to resonance with the leads can 
be designed to produce an effective DoS having sharp res- 
onances and/or pseudogaps near the Fermi energy. These 
features are shown to strongly modify the Kondo state, 
resulting in a wide range of DQD behavior, which we 
explore using numerical renormalization-group methods. 

We show that when dot 1 is coupled to the leads only 
through dot 2, the Kondo resonance on dot 1 develops a 
sizable splitting. Unlike magnetic fields, which produce 
similar splittings, the "band filtering" introduced by the 
connecting dot preserves the Kondo singlet ground state, 
and results in a finite Kondo temperature for complete 



screening of the magnetic moment on dot 1. 

A second configuration, involving coherent dot-dot 
coupling via the leads, can mimic an Anderson impu- 
rity in a pseudogapped host. The device can be tuned 
by varying gate voltages to a QCP separating Kondo- 
screened and free-local-moment ground states. Such 
DQDs offer an attractive experimental setting for sys- 
tematic study of boundary quantum phase transitions. 

The DQD consists of dots 1 and 2 connected to left 
(L) and right (R) leads as well as to each other, as shown 
schematically in Fig. ^ We focus on situations in which 
dot 1 is tuned to have an odd number of electrons in a 
Coulomb blockade valley, so that it has an unpaired spin, 
while near-resonant transport through dot 2 is dominated 
by a single level [l4| having a dot-lead coupling greater 
than the charging energy, so that the dot can be consid- 
ered noninteracting. Our Hamiltonian is therefore 



+ H £ kc] k(T C 3k(T + [ V V a \a C 3ko + H.C.J , (1) 

where a\ a creates a spin-cr electron in dot i (= 1,2), 
ni a = aj^aiu, and cj k(T creates a spin-cr electron of 
wavevector k and energy e k in lead j (= L,R). For 
simplicity, we take the dot-lead couplings Vij to be k- 
independent. We further assume Vm = Viz,, in which 
case the dots couple to the leads only in the symmetric 
combination c^ a = (c^kcr + c^ker) / 'V% and Eq. de- 
scribes double dots effectively coupled to a single lead, 
with Vi = \/2ViL- Near-symmetric couplings can be 
achieved experimentally by appropriate tuning of the 
dot-lead tunneling gate voltages (see, e.g., 0]). Vi and 
the dot-dot coupling A are taken to be real and positive. 

The Green's function (GF) for dot 1 is Gn(w) = {{a lcr : 
a\ a )) = (1 + UiT 11{uj) )G { ^{uj), where u is the energy 
relative to the common chemical potential fi = in the 
leads, rn(w) = {{ni^ a aic, : a\ a )), and G^{uj) is the 
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FIG. 1: (color online) Schematic of the DQD system. Dot 1 is 
Kondo-like (— ei,ei + Ui ^> Ai = TvpoVi , where po is the lead 
DoS), while dot 2 can be treated as a single, noninteracting 
(U2 <C A2 = irpoVi) resonant level. 



noninteracting GF for dot 1 in the presence of dot 2: 



. (2) 



= [uj — Si — ^ 2 /( w ~ £ k)] 1 being the nonin- 
teracting GF for dot i in the absence of the other dot. 

Hereafter, we assume a constant DoS po in the leads. 
In the wide-band limit (half bandwidth D 3> |w|), we can 
formally write [Gn(w)]~ = uj — e\ — + A(uj) + 

iA(u>), where S* 1 (o;) is the proper self-energy, A(u>) = 
tt P2 (uj) [(A 2 - AiA 2 )(w - e 2 )/A 2 - 2\y/A7A 2 ~\ , and 



A(w) = TrpaH A+(w-e 2 )VAi/A 2 



(3) 



with Ai = TrpoVf and p 2 M = A 2 /{tt[(w - e 2 ) 2 + A^]}. 

All information on the coupling of dot 1 to the leads 
and to dot 2 enters Gn through A(lo) (which essentially 
renormalizes the single-particle energy £1) and, more im- 
portantly, through A(uj). This provides a mapping of 
the DQD onto a single Anderson impurity coupled to 
a Fermi system with an effective hybridization function 
A(uj), which "filters" the band states seen by dot 1 and 
modifies its coupling to the leads. We have solved this 
interacting model using an extension of the numerical 
renormalization-group (NRG) method [lal designed to 
handle arbitrary conduction band shapes |fj. 

In order to understand the effects of the nonconstant 
effective hybridization A(w), we consider (i) a side-dot 
configuration, in which dot 1 is coupled to the leads only 
through dot 2 (A 7^ 0, Vi = 0); (ii) a parallel config- 
uration, in which interdot interactions take place only 
indirectly via the leads (A = 0, V\ 7^ 0); and (iii) a more 
general fully connected configuration (A ^ 0, V\ ^ 0). 
(i) In the side-dot configuration, the effective hybridiza- 
tion A(a;) = irp 2 (uj)\ 2 , so the system maps onto an An- 
derson impurity coupled to a Lorentzian DoS. The case 
£2 = 0, which places the peak in A(uj) at the Fermi en- 
ergy uj = 0, yields the most striking properties. 
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FIG. 2: (color online) Side-dot configuration, (a)-(c) Dot-1 
spectral density for Ui = —2e\ = 0.5, £2 = 0, and A2 = 0.02. 
Splitting appears in 4u(w) for larger A such that Tk > 
A2/V2. (d) Effective moment pi vs T. (e) T K vs A. For 
large A, Tk oc A (dashed). (Energies and Tk in units of D.) 



Figure [3 presents results for e 2 = and U\ = — 2ei, 
for which parameters the model exhibits strict particle- 
hole (p-h) symmetry. Figs. Ela)-(c) show the spec- 
tral density 4 u (w) = -Im Gn for A 2 = 0.02D 
and different interdot couplings A. For small A [Fig. 
EJa)], the spectral density resembles that for a constant 
DoS, its main features being broad Hubbard bands cen- 
tered near uj — E\ and s\ + Ui, and a sharp reso- 
nance at uj — having a width of order the Kondo 
temperature Tk (defined below). In this regime, the 
nonconstant A(w) manifests itself through a generalized 
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yj/[7rA(0)], where 



Fermi-liquid relation 

ip = /"^ [(dA/do;)ReGii(w) - (dA/<L))ka.G n (u)]du), 
i.e., An(0) is smaller by a factor of cos 2 <p than the stan- 
dard result T] for a flat band with the same A(0). 

For larger A, such that Tk > A2/V2, the spectral 
density is qualitatively different. The Kondo resonance 
initially rises under the influence of the relatively weak 
hybridization found for \uj\ > Tk- However, the upturn 
in A(uj) at \uj\ < A 2 causes Ah(uj) to drop to satisfy 
-^11 (0) < l/[7rA(0)] (see above), splitting the resonance 
into two distinct peaks [Fig. Hfb)]. As A increases fur- 
ther, the dip deepens and the Kondo peaks move out, 
eventually subsuming the Hubbard bands [Fig. Etc)] . 

The splitting of the Kondo peak and the suppression 
of .Aii(O) might be supposed to signal the destruction 
of the Kondo singlet ground state (as is the case in a 
magnetic field). However, this interpretation is refuted 
by the NRG many-body spectra, which show that the 
Kondo ground state is reached even for Tk 3> A 2 . 

The progressive screening of the localized spin with 
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decreasing temperature T can be seen Fig. Hid), which 
plots the square of the effective free moment on dot 1, 
p\(T) = Txj mp , where Xi mp is the dot contribution to 
the zero-field susceptibility. The behavior for Tk < A 2 
follows that for a constant DoS 0, crossing at T w 
from the free-orbital regime (/zf ps 1/8) to the local- 
moment (LM) regime (p.\ ps 1/4), and then for T -C 
to the strong-coupling (SC) Kondo limit in which the 
magnetic moment is totally screened (/xi = 0). 

For Tk > A 2 , the system still reaches the Kondo fixed 
point with jtt| = as T — ► 0. However, the dot spin now 
exhibits an interesting window of diamagnetic behavior 
(Ximp < 0), which becomes more pronounced as Tk in- 
creases [Fig. EJd)]. This diamagnetic region arises from 
a negative term cx dp2(co)/dw in Ximp |l6| . 

Since /i 2 in all cases passes from 1/8 at high temper- 
atures to at T = 0, we define the Kondo temperature 
using the standard criterion (4(T K ) = 0.0701 [ljj. T K , 
shown in Fig. |2Ie) for three different A 2 values, increases 
rapidly for small A, and satisfies T% oc A in the noninter- 
acting narrow-band limit A 3> A 2 , Ui/2. 

The splitting observed in the Kondo resonance and the 
diamagnetic region in p\{T) are produced by the sharp 
peak in p 2 (w) at u> = resulting from the resonance 
in dot 2. The resonance acts as a filter for the higher- 
energy states in the leads, reducing the effective conduc- 
tion bandwidth connected to dot 1. This interpretation, 
which is consistent with similar findings of a negative 
Ximp in narrow-band systems clearly differentiates 
the side-dot behavior from the peak splittings due to co- 
herent coupling to a second Anderson impurity and the 
resulting suppression of the singlet state [llj . 

The side-dot behavior can also be understood as arising 
from interference between resonances. This can be seen 
by considering the noninteracting spectral density = 
-hnG$/n for e x = £ 2 . For A < A 2 /V2, A$ has a sin- 
gle peak (width ~ A 2 ) at ui = e 2 ; whereas for A > A 2 /v / 2, 
there are two peaks at u> = e 2 ± yj\ 2 — A 2 /2, arising 
from interference between the ui — single-particle reso- 
nances on the two dots. The NRG results for the inter- 
acting case are closely analogous: For Tk < A 2 /\/2, An 
has a single peak (width ~ Tk) at ui — £ 2 = 0, while 
for Tk > A 2 /a/2, there are two peaks, in this case re- 
sulting from interference between the u> — many-body 
Kondo resonance and the u) = single-particle resonance 
in p2- The separation of the peaks in An increases with 
Tk — A 2 , further supporting the analogy, 
(ii) The parallel configuration (A = 0) exhibits very dif- 
ferent behavior. For £ 2 = 0, the hybridization A(w) [Eq. 
0] vanishes at u> — as \uj\ t with r — 2 [see Fig.^Ja)], 
and the DQD setup maps onto an Anderson impurity in 
a pseudogapped host . The properties of such system 
depend strongly on the exponent r and on the presence 
or absence of p-h symmetry 0, Q. For r — 2, the SC 
phase is inaccessible at p-h symmetry, but away from 
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FIG. 3: (color online) Parallel DQD with Ui = 0.5, Ai = 0.05, 
A2 = 0.02, and £2 = 0. (a) Hybridization A(w) vanishes as 
|u;| 2 at u — 0. (b)-(d) Dot-1 spectral density Aii(u), and 
(e) effective moment /nf(T), for various ei. For ei ~ —Ui/2, 
no Kondo effect occurs: Au(0) = in (b), and /x?(0) = 1/4 
[■ and □ in (e)]. This LM phase is separated by a QCP at 
£1 = e* > —Ui/2, characterized by a featureless Au(w) in (c) 
and (ii(Q) = 1/6 [• in (e)], from the SC phase [(d), a and a in 
(e)] in which An vanishes at ui = and peaks at some u> > 0. 



this special limit, a QCP separates SC and LM phases. 

In the p-/i-symmetric case E\ = —Ui/2, the Kondo 
resonance in An disappears completely [Fig. |3fb)], but 
Hubbard bands are still present (arrows in inset). An 
vanishes at w = as uj 2 , and = 1/4 [Fig. E£ e )]j 

as expected in the LM phase |a] where no Kondo effect 
occurs. When p-h symmetry is broken by increasing £1, 
the same qualitative behavior persists [squares in Fig. 
'6{c)\ until the QCP is reached at £1 = e*. where An is 
nearly featureless [Fig. E{c)] and /u 2 (0) = 1/6 [circles in 
Fig. 03[e)]. For S\ > £*, the system enters the SC phase 
in which /i 2 (0) = [triangles in Fig. Efe)], indicating 
complete Kondo screening of dot 1, and An again goes 
to zero as ui 2 at u> — 0, but (in contrast to the LM phase), 
there is a distinct peak at positive u> [Fig.^d)]. 
(iii) In the fully connected configuration (V\, V2, and A 
all nonzero), A(w) has an asymmetric Fano-like shape, 
peaking at £ 2 + y / A]~Af/A, and vanishing as (lo — loq) 2 at 
lu = £ 2 — A-\/A 2 /Ai. Here, the DQD properties can be 
controlled not only by tuning A, Ai, or A 2 [as in (i) and 
(ii)], but also by using external gate voltages to vary £ 2 . 

Figures|H[a)-(e) show the evolution of A(w) as £ 2 shifts 
across the Fermi energy. Generically, the dot-1 spec- 
tral density features a p-/i-asymmetric Kondo resonance, 
whose width Tk is primarily determined by A(0) [Figs. 
EJf),(h),(j)]. However, when the peak of A(w) approaches 
oj = 0, Tk is enhanced and the Kondo peak splits [Fig. 
Etg)]. By contrast, when the zero of A(ui) approaches 
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FIG. 4: Fully connected DQD with Ui = -2si = 0.5, 
Ai = A 2 = 0.02, and A = 0.1. (a)-(e) Hybridization A(w) for 
various £2, and (f)-(j) the corresponding Aii(w). An asym- 
metric Kondo peak appears when A(w) is featureless near 
u = (a, e). The peak splits as the resonance in A(lu) ap- 
proaches w = (b), reappears in the intermediate region (c), 
and disappears altogether when a pseudogap forms (d). 



the Fermi energy, the Kondo temperature decreases, and 
if ujq is tuned to zero, pseudogap behavior is recovered 
with complete suppression of the Kondo peak [Fig.^i)]. 

In summary, we have shown that double quantum-dot 
(DQD) systems with one of the dots in the Kondo regime 
can be tailored experimentally to explore the effects of a 
nonconstant density of states (DoS) on the many-body 
ground-state properties. In setups where the Kondo dot 
is decoupled from the reservoirs, the Kondo resonance 
on that dot undergoes zero-field splitting for large inter- 
dot coupling. This can be understood as the coherent 
interaction between the many-body Kondo state and a 
single-particle resonance in the second dot. Although 
the Kondo peak in the spectral density at the Fermi en- 
ergy is suppressed, the Kondo singlet state is robust, and 
the localized spin is completely screened at low temper- 
atures. In this regime, the system also passes through a 
temperature window of diamagnetic behavior, similar to 
that seen in narrow-band systems. 

For weak interdot couplings, the Kondo state is sup- 
pressed by the presence of a pseudogap in the effective 
DoS, and a quantum phase transition takes place between 
local-moment and Kondo-screened phases. A quantum 
critical point on the boundary between these phases can 
be reached by appropriate tuning of the experimental 
parameters. Thus, DQD systems provide a rare example 
of a controlled setting in which to investigate quantum 
critical behavior in a strongly correlated system. 
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